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Previous work on optimum gain control is extended to an important 
class of diversity receivers used for digital data transmission through jading 
media and for radar. As in the single diversity case the optimum gain 
{which yields minimum average cost of receiver saturation) is extremely 
insensitive to relative costs of saturation at the upper and lower dynamic 
range bounds. The sensitivity to relative cost decreases as the order of 
diversity increases. 

Optimum gain and performance characteristics are given from which 
dynamic range requirements for diversity receivers can be deduced. 

I. INTRODUCTION 

A good part of detection theory literature deals with the determina- 
tion of statistically optimum or near optimum receiver structures. 
However, in any practical implementation of these receivers the signal 
processing must be performed by components of finite dynamic range. 
To effectively use the amplitude range of a signal processing chain 
it is common to scale the received signal by adjusting the receiver 
gain. Optimum gain settings for minimum average cost of excluding 
(from a receiver's finite dynamic range) the envelope of a narrowband 
signal plus gaussian noise were presented last year. 1 Here similar 
results are presented for an important class of diversity receivers 
used for comunications through fading media and for radar. For 
the single diversity {M = 1) case, these results reduce to those 
given previously. 

II. PREVIOUS RESULTS 

Consider the problem of determining the normalized attenuation, 
a to optimally scale a positive homogeneous functional, £, of the 
received signal so that the average cost, I, of excluding £ from the 
receiver's dynamic range is minimized. It follows from Ref. 1 that 
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the average exclusion cost is given by 

i= [' p.® <% + » fV©# (i) 

Jo Jnd 

in which d denotes the dynamic range of the receiver (such that D(dB) = 
20 logio d, d > 1), v is the ratio of cost of saturation at the upper dy- 
namic range bound to the cost of saturation at the lower, co is a vector 
parameter determined by signal noise and channel conditions, and p u 
is the probability density function of £. When the optimizing value of a 
is a stationary point of I it can be found as a real positive solution to 

p a (a) = udp w (ad). (2) 

If a 8 is the optimum a then the minimum average exclusion cost is 

I = p u ( a ,) + v [l - P u (a, d)] (3) 

in which P u is the cumulative distribution corresponding to p a . For 
v — 1, I becomes the exclusion probability. Ref. 1 considered this prob- 
lem in detail for the case in which £ represents the envelope of a narrow- 
band signal plus gaussian noise received through a Rician fading me- 
dium. The results are based upon the solution of (2) for the case in 
which p a is the Rician 2 probability density function defined by 

Pt(0 = 6 exp [- (f + t 2 )/2]/o(t£) 7, k^ (4) 

where y is a suitably defined signal-to-noise ratio. (y<iB = 20 logioy). 

III. GAIN SETTINGS FOR DIVERSITY RECEIVERS 

In various diversity receivers formation of the test statistic leads 
to the generalized Rician probability density function given by 

p.(JZ) = R(R(M)*/y)"- 1 exp [-(R 2 + -? /M)/Wm-xW{M)*\ (5) 

R,y, ^0 
M = 1, 2, • • • 

where I K denotes the modified Bessel function of the first kind and order 
K and w is the vector (7, M ). Such is the case for example in square-law 
combining ilf-fold diversity receivers for noncoherent frequency shift 
keyed signaling through Rician or Rayleigh (if 7 = 0) fading channels, 
in radar receivers using post detection square law integration of M 
pulses, f and in partially coherent diversity reception of iV-ary orthogonal 

t In these cases the functional £ is the test statistic. More generally however, 
the functional used for determining the optimum gain need not be actually 
formed in the receiver. 
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signals transmitted through M independent slow Rician fading chan- 
nels. 3,4 
The probability density function (5) has interesting properties. It 
can be shown that (5) is the probability density function of the 
square root of the sum of squares of M independent normalized 
Rician variates, each variate having a probability density function 
of the form (4) with y replaced by y/(M)%. For M = 1, (5) of 
course reduces to the Rician probability density function (4), and 
if in addition y = it becomes the Rayleigh probability density 
function. The density (5) can be viewed as the probability density 
function of the square root of a noncentral chi-square variate with 
2M degrees of freedom and noncentrality parameter y 2 . With y = 
it becomes the density of the square root of a chi-square variate with 
2M degrees of freedom. 

In many practical cases y 2 is proportional to the ratio of the total 
specular energy received via the M diversity branches to the sum of 
the scatter and noise energy received via any diversity branch (as- 
suming that this latter sum is the same for any diversity branch). 
Thus y 2 /M can be thought of as the power signal-to-noise ratio per 
diversity branch or per pulse in the case of time diversity if, as is 
commonly assumed, the diversity branches are statistically independ- 
ent but have identical parameters. 

Since (5) arises in various applications, let us consider the canonical 
problem. Specifically the solutions to (2) will be obtained in which p u 
is given by (5). Letting o = A(2)* and a = y/(M) { leads to the following 
transcendental equation for A : 

A 2 = A;+ [(M + l)/2] Al + [l/(d 2 - 1)] 

• {hi / JM M d(2) i ] - In 7jr-iM(2)*]} 

(6) 



in which 



2 A 2\rxd (7) 

A - d2 _ 1 KO 

determines the optimum required attenuation for the Rayleigh case 
with unity cost ratio (v = 1) , and 

Al ^ J^r- (8) 

a — 1 

For the single diversity case (M = 1), (6) reduces to the trans- 
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cendental equation encountered previously. 1 One is thus led to seek 
an iterative solution along the same lines. Before obtaining the 
required iteration equations consider some properties of (6). The 
equation can be written in the form 

A 2 = A; + MAI + [l/(d 2 - l)]{ln [M d(2) i y t ~ 1 hr-i(aA d(2) i )) 

- (M - 1) In d 2 - In [(aA(2)Y~ 1 Iu-^A(2) i ))} . (9) 

Combining the logarithmic terms on the right side of (9) and using 
the fact that I K (Z) -» (Z K )/2 K K\ as Z -> it is seen that the quan- 
tity in braces goes to zero as y -» 0. Hence the optimum attenuation 
for the chi-square case (y = 0) is determined explicitly by A = A cs 
where 

Al 4 A 2 + MAI . (10) 

In the same manner it is seen that if A 2 C , = 0, then A = is a solution 
to (9) for any y and d. It is easy to show that the right hand side of (9) 
is an even function of A having a minimum of A 2 , at A = 0. The left 
hand side of (9) is of course a standard parabola centered at the origin. 
These curves (i) do not intersect if A 2 C , < 0, (it) intersect only at A = 
if A 2 C , = 0, and (in) intersect at positive (and negative) values of A if 
A\, > 0. Thus meaningful values of A which minimize I are stationary 
if and only if A 2 C3 > 0. From (7), (8), and (10) this requires vd 2M > 1, 
which for M = 1 reduces to the constraint encountered previously. 1 
The solution to (6) or (9) can be obtained using the extrapolated 
iteration scheme described in Ref. 1. The iteration formulas require the 
derivative of the right side of (9) which can be found using the identity 

d/dt [rUm = rV-i(fl n - ■ ■ ■ -2, -1, 0, 1, 2, ■ ■ • . (11) 

The result is 

WA* y(2)* I I M - 2 [Ada(2)t) I M - t [Aa<&]\ (m 

n* ) - 2(d 2 _ 1)A \a Im _ i[A da(2)i] Ilf _ l[Aam j w 

in which / denotes the right hand side of (9) and the prime denotes dif- 
ferentiation with respect to the argument. 

For computational purposes it is convenient to define the functions 

¥xfr) A (exp -f) J*(f) (13) 

which are uniformly bounded on the semi-infinite interval [0, oo]. 
For any argument £ these functions can be readily generated by 
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numerical techniques using the recurrence relations and asymptotic 
expansions for the modified Bessel functions. 

Using (12) and (13) the following iteration formulas to solve 
(6) are obtained. 

A 2 i+1 = F t -fc(Aj - F<) ft* -1 

F t = A: + [(M + l)/2] Al 



+ fff + M* - 1)] 



' IjMJij <fa(2)») 



(14) 



a(2)» / , *„- 2 rA,-rfa(2)*] _ ^- 2 [.4,-a(2)^ 

*" " 2(d 2 - 1)A, T *jr-iU, da(2)*} ^,-,[i,a(2) ! : 

0. = G,/(l - Gi) Gi ?* 1 

The iteration is begun with i = 1, 7 small, and A\ = A;;, and stopped 
when \ (A i+1 — Ai)/Ai \ is less than the allowable error. By this method 
the optimum required normalized attenuation was found for various 
values of v, 7, d, and M. 

Inasmuch as the optimum attenuation satisfies the nonlinear equation 
(6) it will be helpful in interpreting results to find some useful approxima- 
tions. Accordingly one notes that for y 2 /M >>> 1 and c? >$> 1 the second 
term in the brackets on the right side of (6) is negligible compared with 
the first. Then taking Im(x) tt exp x leads to a quadratic equation in A 
whose solution 

approximates the required attenuation over the range specified. 

Fory*/M « 1, on the other hand, one may take I M {x) W (x/2) M /M\ 
in (6). Some further approximations and manipulation lead to the 
result 

A « A r .[l + (y 2 /2M)] h r/M « 1 (16) 

which is exact for 7 = 0. 

When the solution to (6) is found for given parameters the mini- 
mum average exclusion cost can be determined. For the probability 
density function (5), (3) becomes 

1=1- Q M [a, .4(2)*] + pQ M [a, Ad(2) i ] (17) 
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where 

Qm(x, y) = P rCf/x)"- 1 exp [-(f 2 + .T 2 )/2]i, w _,(xf) rff (18) 
is the generalized Marcum ©-function. 8 

IV. EFFECT OF COST RATIO ON REQUIRED ATTENUATION 

Since one cannot easily decide how much better or worse it is for 
the receiver to saturate at its upper limit than at its lower limit, the 
ratio, v, is generally difficult to assess accurately. It is interesting (even 
fortunate!) that here as in the single diversity case, 1 the solutions 
obtained using (14) show that the optimum receiver attenuations 
for a wide range of cost ratios do not differ appreciably from those 
for the minimum exclusion probability case {v = 1) . 

For given v, y, d and M one may define the sensitivity, S , of the 
optimum attenuation to the cost ratio by the difference in required 
attenuation between the given case and the corresponding minimum 
exclusion probability case. Specifically, the sensitivity to cost ratio is 

S c (p, 7, d, M) 4 20 log 10 A{v, y, d, M) 

- 20 1og ]0 ^(v = 1,7, d,M) (19) 

in which the functional dependence is shown explicitly. In (19) a positive 
value of S c indicates an increase in required attenuation compared with 
the minimum exclusion probability case. It can be shown that (i) the 
sign of (19) depends only on v (positive if v > 1, negative if v < 1), and, 
iii) | S e (v, 7, d, M) \ ^ | S c (v, 0, d, M) \ . Thus, one can define a maximum 
sensitivity 

S* 4 S e (v, 0, d, M) = 10 log 10 A], - 10 log 10 MAl (20) 

where the second equality follows from (19) and (10). Using (7) 
and (8) , (20) can be written 

S* = 10 log 10 [1 + (W2MD)] (21) 

in which 

v dB ^ 20 log 10 v (22) 

is the cost ratio expressed in dB. S* c is an easily calculated bound which 
gives, with the correct algebraic sign, the maximum change (in dB) of the 
optimum required attenuation from the optimum for the minimum ex- 
clusion probability case. Figure 1 is a plot of S*. It follows from (21) 
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Fig. 1 — Maximum change in required receiver attenuation caused by nonunity 

cost ratio. 



that for cost ratios in the range —MD ^ v dB ^ 2MD, the maximum 
change in required receiver attenuation is less than about 3 dB. Equiva- 
lently, with y, v, and d fixed, the maximum sensitivity, aS* decreases as 
the order of diversity, M, increases. 

Since the optimum attenuation is extremely insensitive to cost 

ratio for typical parameters, the minimum exclusion probability 
case (v =1) is of special import among all average cost criteria of 
the form (1). The numerical results presented in this paper, there- 
fore, include only the case v = 1, although the formulas derived apply 
more generally and can be used to generate numerical results in an 
entirely similar manner. 



V. EFFECT OF DIVERSITY ON REQUIRED ATTENUATION 

It is interesting to consider how the order of diversity affects the 
optimum required attenuation. Accordingly, in a manner analogous 
to (19) one can define the difference in required receiver attenuation 
resulting from diversity by 

S m (v, y, d, M) = 20 log I0 A(v, y, d, M) 

- 20 log 10 A(p, 7, d, M - 1). (23) 
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Let 

A n (v, y, d, M) 4 20 log 10 A(p, t, d, M) (24) 

be the required normalized attenuation in dB. Then (23) can be 
written 

S m (v, 7, d, M) = A n (v, 7, d, M) - A n (v, 7, d, 1). (25) 

Values of <S,„(1, y, d, M) were obtained for various y, d, and M 
using (14) and (23). These are shown in Fig. 2 where all quantities 
except M are in dB. The optimum normalized attenuation A„(l, y, 
d, M) required for the minimum exclusion probability case can be 
found using (25). Specifically one finds S OT (1, y, d, M) from Fig. 2, 
and adds to it the quantity A„(l, y, d, 1) from Fig. 6 in Ref. 1. 
Notice that in Ref. 1 only the single diversity case (M = 1) was 
considered so that the functional dependence of A„ on M was sup- 
pressed in the notation. That is, A„(v, y, d, 1) here is identical to 
A„(i', y, d) in Ref. 1. 

From Fig. 2 it can be seen that if y is sufficiently small (or large) , 
S m is positive (or negative) so that more (or less) attenuation is 
required if multiple diversity is used than would be required if the 
same specular energy were concentrated in a single diversity branch 
or pulse. Also for sufficiently small (or large) y the required attenua- 
tion increases (or decreases) as the order of diversity, M, increases. 
There is of course a transition region which bridges the above cases 
and in which, for y fixed, the differences S m cross one another depend- 
ing on the particular values of M, and D (and, in the general case, 
v). The curves for y dB = 15, for example, exhibit this behavior. 

Using (10) and (16) it can be shown that for y -» 

S m (v, y, d, M) « 10 log 10 M 

■ in l«. f l + fr dB /2MZ)) l ■ inw [ l + fr721f) 1 (m 
+ 10 l0gl ° L 1 + (»Mb/22» . + 10 l0g, ° L 1 + (t72) J (26) 

which is exact for y = 0. For the minimum exclusion probability 
case and y = (26) yields S m (1, 0, d, M) = 10 log 10 M which is in- 
dependent of d. Similarly it can be shown using (15) that 

lim S m (", y, d, M) = -10 log 10 M (27) 

which is independent of v and d. The differences S m for y dB - =fc °° 
therefore appear as horizontal fines in Fig. 2. One also observes that 
over the range of parameters shown, the limit (27) is approached within 
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Fig. 2 — Differences in required optimum attenuation 
resulting from diversity, v = 1. 
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0.5 dB for y dB = 30. Noting that for y° < 7 < 7*, S m is bounded by 

S m (v, 7 * 5 d, M) < S m (v, 7, d, M) < S m (u, y°, d, M) (28) 

it follows from (26) through (28) that for the minimum exclusion proba- 
bility case all the differences S m (l, 7, d, M) lie between two horizontal 
lines in Fig. 2 determined only by the order of diversity. | S m (l, 7, d, M) | ^ 
10 log 10 M. The difference between the required optimum attenuation 
for dual diversity and that required for single diversity with the same 
total received specular energy is less than about 3 dB. 

VI. EXCLUSION COSTS FOR DIVERSITY RECEIVERS 

The optimum normalized attenuations obtained using the iteration 
equations (14) were used to obtain the minimum average ex- 
clusion costs (17) for the case v = 1. These are shown in Fig. 3(a) 
and for smaller values of D in Fig. 3(b). The generalized Q function 
(18) was evaluated by computer, using relations derived from those 
given by Sagon. 5 

It can be seen that for small values of y, the smallest dynamic range 
D required to obtain a given exclusion probability decreases rapidly 
as the order of diversity is increased; the most substantial decrease is 
obtained in going from single to dual diversity. This trend is lessened 
as the available signal-to-noise ratio y increases. As a matter of fact 
if y is sufficiently large (for example, y = 20 dB) the dynamic range 
required to achieve a given exclusion probability increases as M 
increases. However at the large values of y where this latter effect is 
apparent, even modest values of D yield extremely small exclusion 
probabilities. Moreover on the types of channels where diversity re- 
ceivers are useful one would generally encounter small values of y. 

Consider a diversity receiver operating in a small signal-to-noise 
ratio and let the dynamic range of the components used be such that 
the probability of excluding the signal at any point in the receiver is 
the same throughout. Then it follows from Fig. 3 and the foregoing 
discussion that the dynamic range required of the components used 
in the post-combining portions of the receiver may be considerably 
smaller than that required of those components used in the individual 
diversity branches. 

VH. SUMMARY AND CONCLUSIONS 

An important class of diversity receivers used for communications 
through fading media and for radar is considered. The required gain 
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Fig. 3 — Minimum exclusion probabilities for diversity receivers. 
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is determined which minimizes the average cost of excluding from a 
finite dynamic range the signal appearing in the post-combining por- 
tions of the receiver. For the single diversity case {M =1) the results 
reduce to those given previously. 

It is shown that the required receiver gain is extremely insensitive 
to the relative costs of saturation at the upper and lower dynamic range 
bounds, differing at most by about 3 dB from the optimum for the (equal 
cost) minimum exclusion probability case for relative costs in the range 
— MD ^ v dB ^ 1MD. One also finds that the sensitivity to relative cost 
therefore decreases as the order of diversity increases. 

The difference between the required optimum receiver gain for 
various orders of diversity M , and that required for a single diversity 
receiver having the same total received specular energy is considered. 
Exact differences are given for the minimum exclusion probability 
case, and it is shown that these are less than 10 logm-M dB independent 
of other parameters. Bounds on the difference are also given for non- 
unity cost ratio. 

Performance characteristics derived show minimum exclusion prob- 
abilities obtainable as a function of dynamic range for various signal- 
to-noise ratios and orders of diversity. For a small signal-to-noise 
ratio the dynamic range required of the components used in the post- 
combining portions of the receiver can be considerably smaller than 
that required of those components in the individual diversity branches 
in order to achieve uniform exclusion probability throughout. 

Notice that in some applications the normalization assumed in 
writing (4) and (5) may depend upon M and y. This fact must be 
accounted for if one is calculating the actual required attenuation 
from the required normalized attenuation discussed in Sections IV 
and V. The optimum exclusion costs however, depend on the normal- 
ized attenuation and not on the normalizing factor. The results of 
Section VI therefore apply directly. 
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